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We study the electroweak phase transition in an effective 3-dimensional theory for a Higgs mass of about 70 
GeV by Monte Carlo simulations. The transition temperature and jumps of order parameters are obtained and 
extrapolated to the continuum using multi-histogram techniques and finite size analysis. 



1. Introduction 

One approach to lattice calculations of the elec- 
troweak transition is based on an effective 3- 
dimensional S'[/(2)-Higgs model. It is attrac- 
tive phenomenologically because it circumvents 
the problem of putting chiral fermions on the lat- 
tice. Due to dimensional reduction, fermions as 
well as non-static bosonic modes contribute to 
the effective action. In contrast to QCD, dimen- 
sional reduction should work for the electroweak 
theory around and above the transition temper- 
ature because g"^ is small. For the electroweak 
phase transition this approach has been pioneered 
by Farakos et al. (see e.g. ||]). In its simplest 
version the dimensionally reduced effective the- 
ory is again an S'C/(2)-IIiggs theory with just one 
doublet. 

2. The lattice model 

On the lattice, we study the SU{2)-Eiggs sys- 
tem with one complex Higgs doublet of variable 
modulus. The gauge field is represented by uni- 
tary 2x2 link matrices Ux,on the Higgs fields are 
written as = PxVx {V^ G SU{2)). The lattice 
action is 
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(pI = iTr Up denotes the SU{2) pla- 

quette matrix), with 



gl Pg 8 ^SOGeVy Pa' 
2(1 - 2/3fl) 
6 -t- a^rrii 



(2) 

(3) 
(4) 



(a is the lattice spacing). The lattice model de- 
fined by (0) is numerically studied at given cou- 
plings Pg,Ph and A3/5I. 

In the search for the phase transition, bulk vari- 
ables like 



link 



j-rlll^' ($^C/.,a<i>.+o) (6) 



LxLyLz 



SLxLyLz 



are used. 

The update is a combination of 3d and -id 
Gaussian heat bath for the gauge and Higgs 
fields, respectively, and Higgs reflections. Most 
of the Monte Carlo data have been obtained on 
QUADRICS parallel computers. 



* Contribution presented by M. Giirtler 
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3. Phase separation, equal weight and 
mixed phase configurations 

As usual, the search for the phase transition 
point requires extensive application of the multi- 
histogram technique 

We have studied the phase transition driven 
by 7713. Then the lattice Higgs self-coupling (3^ 
varies with I3h (see (Q)). Therefore, the reweight- 
ing uses not only Eunk , but and p'^ at the same 
time. 

We have determined the finite volume pseudo- 
critical [3hc{L) by the minima of the Binder cu- 
mulants and the maxima of the susceptibilities 

B^.{L,Ph) = 1-^7^^, (7) 

C,.{L,Ph) = {{p'f)-{p'f (8) 

(in the same way for other observables) as well as 
using the equal weight method. 

The application of the equal weight criterion 
requires a procedure to separate the (measured 
or reweighted) histogram into contributions at- 
tributed to the pure phases. In addition, there 
are inhomogeneous (mixed) configurations con- 
tributing to the histogram. Our main assumption 
is that the pure phases can be described by Gaus- 
sian distributions for any volume averaged quan- 
tity. The normalised histogram has been repre- 
sented as a weighted sum of three histograms 

P{P^,Ph) = WbPb{p^,PH) + WsPs{p^,fiH) 

{p\Ph) (9) 

with Wf, -\-Ws +Wmix = 1- Wb,s dcnotcs the weight 
of the broken/symmetric phase, Wmix that of the 
mixed state. 

The positions, widths and weights of the pure 
phase histograms at a given have been ob- 
tained by fitting the outer flanks of the two- 
peak histogram to Gaussian shape. This fixes 
the weight Wmix and the p^ distribution to be at- 
tributed to configurations with domains of both 
phases in equilibrium. The pseudo-critical Phc is 
found according to the requirement Wb = Wg. 

Another kind of phase separation was used to 
estimate the jump of the plaquette and correla- 
tion lengths of "pure" phases at the critical point. 



The aim is to remove successful tunneling escapes 
and unsuccessful tunneling attempts towards the 
"wrong" phase from what should then be consid- 
ered as the Monte Carlo trajectory restricted to 
the "right" phase. 

The procedure rescans the records of p^ which 
has a well separated two-peak signal for all con- 
sidered volumes. Referring to this variable a 
lower cut for the upper (broken) phase and an 
upper cut for the lower (symmetric) phase can be 
chosen. The cuts are determined in such a way 
that the remaining histograms (for the "pure" 
phases) are almost symmetric around their max- 
ima. If the Monte Carlo history of enters the 
range of a certain phase and stays there for a 
number of iterations (larger than the autocorre- 
lation time without tunneling but smaller than 
that with tunneling), the sequence of configura- 
tions is considered to belong to that phase until 
the trajectory leaves it. 

4. Phase Transition Localisation 

To demonstrate the two-peak structure we 
show in Fig. |] the measured histogram of on 
lattices 48^ and 64^, all at /?g = 12, for Ph values 
nearest to the respective pseudo-critical [3hc{L). 
The positions of the maxima change already only 
slightly with the volume. 
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Figure 1. Measured histograms of p^ for /?g = 12 

In Fig. ^ results of the multi-histogram inter- 
polation of our data for Pa = 12 for the Binder 
cumulant of Eunk are presented. Finiteness and 
shrinking of the Binder cumulant with increasing 
volume present evidence for the first order nature 
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of the transition at Higgs mass M'^ = 70 GeV. 
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Figure 2. Binder cumulant for Eunk, Pg = 12 

An equal- weight histogram {f3Hc{L) = 
0.3435441) for a lattice size 64"^ is presented 
in Fig. H together with the Gaussians describing 
the pure phases. The distribution attributed 
to mixed configurations with domains of both 
phases in equilibrium is well identified between 
the two peaks. 




3.2 3.4 3.6 3.8 4.0 4.2 



Figure 3. Equal weight histogram 

The various pseudo-critical (3hc{L) values for 
the different methods applied to Eii„k and are 
collected in Fig. ^, plotted versus for Pg = 
12. Corresponding to each method, a fit 
has been used to yield a respective /3^j,. The 
extrapolations nicely coincide as expected. 

In Table 1 the extrapolations for each method 
are collected together with the average /3^^ for 
(3g = 12 and 16. The result is translated into a 
physical temperature and an "exact" Higgs mass 
Mh using the correspondence to quantities in the 
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Figure 4. Inhnite volume extrapolation of 

Phc(Pg = 12 j 

4-dimensional theory in l~loop accuracy (no fi^- 
dependence, see [||). 

These numbers are given for the case of the 
S'C/( 2 )-Higgs theory without and with fermions 
(i.e. including the top quark with a mass of 175 
GeV). Comparing the temperatures, there seems 
to be not much space left for 0(a) corrections. 
The "exact" Higgs mass is practically the same. 



Pg = 12 


B C Wb = Ws 


Elink 


0.3435434 0.3435430 

0.3435429 0.3435429 0.3435441 




0.3435433(6) 




no fermions mt = 175GeV 


Tc/GeV 
Mh /GeY 


150.94(1) 107.05(1) 
64.77 69.42 




Pg - 16 


B C 


Elink 


0.3407950 0.3407937 
0.3407943 0.3407939 


I3hc 


0.3407942(6) 




no fermions mt=175GeV 


Tc/GeV 
Mn/GeY 


151.27(1) 107.17(1) 
64.77 69.46 


Table 1. Inhnite volume limit for Phc at 
M^j = 70 GeV 

For comparison, at the smaller coupling (M|f = 
35 GeV) the transition temperature (without top) 



is Tc = 76.2(1) GeV with the Higgs mass Mh = 
29.5 GeV. This has been obtained for gauge cou- 
plings in the range from /3g = 12 to 20 on lattices 
of size 40^ and 20^ [|. 
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5. Condensate discontinuities 

The jumps in (p^) and (p*) are connected to 
the renormalization group invariant discontinu- 
ities of the quadratic and quartic Higgs conden- 
sates. The two-state signal for (p^) and {p^) 
is still clearly visible for all lattice sizes consid- 
ered at the higher Higgs mass of = 70 GeV, 
where the transition turns out much weaker than 



at M 



H 



35 GeV. 



Corresponding to the different criteria applied 
for the definition of the pseudo-critical Phc{L) we 
obtain histograms of the various operators just at 
the respective pseudo-criticality. A collection of 
discontinuities of (p^) for various finite lattices 
(read off from the maxima of the corresponding 
histograms) is shown in Fig. |^. As in the analy- 
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Figure 5. InEnite volume extrapolation of A(p^) 
and A(p4) 

sis of Ref. Q we found that the size dependence 
of the jumps for all available lattice sizes is best 
described by a 1/L^ fit. The extrapolation to in- 
finite volume is given in Table 2 in lattice units. 





Bp2 Cp2 average 


/3g = 16 


0.490(9) 0.479(8) 0.485(6) 
0.372(8) 0.362(7) 0.367(5) 




Bpi Cpi average 


/3g = 12 
I3g = 16 


4.86(9) 4.81(8) 4.84(6) 
3.62(8) 3.51(7) 3.56(5) 



Table 2. Infinite volume limit for A{p^) (upper) 
and A{p*) (lower part), = 70 GeV 

Additionally, the expectation value of the aver- 
age plaquette (P) shows a discontinuity as well. 
This jump is numerically a tiny effect at the larger 



Higgs mass. Nevertheless, we are able to estimate 
it using the phase separation technique discussed 
earlier. 

In Table 3 the jump A(P) is reported for 
values nearest to the critical ones at = 70 and 
35 GeV at lattice sizes 64"^ and 40'^, respectively. 





A(P) 


= 70 GeV (3g = 12 
/3g-16 


0.00037 
0.00015 


= 35 GeV /3g = 12 


0.00370 



Table 3. Estimated plaquette jump A{P) 

The relation of the measured quantities to con- 
tinuum physics as well as further consequences 
are discussed in the related contribution of A. 
Schiller §. 

6. Summary 

We have studied the electroweak phase tran- 
sition for Higgs masses up to 70 GeV. Its first 
order nature has been demonstrated. The criti- 
cal parameters were determined using the multi- 
histogram technique in conjunction with the 
equal weight and other more standard criteria 
giving consistent results. The jumps in several 
quantities are reported and extrapolated to the 
continuum. 
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